Programmable unknown quantum-state discriminators with multiple copies of 
program and data: A Jordan basis approach 
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The discrimination of any pair of unknown quantum states is performed by devices processing 
three parts of inputs: copies of the pair of unknown states we want to discriminate are respectively 
stored in two program systems and copies of data, which is guaranteed to be one of the unknown 
states, in a third system. We study the efficiency of such programmable devices with the inputs 
prepared with n and m copies of unknown qubits used as programs and data, respectively. By 
finding a symmetry in the average inputs, we apply the Jordan basis method to derive their optimal 
unambiguous discrimination and the minimum-error discrimination schemes. The dependence of 
the optimal solutions on the a prior probabilities of the mean input states is also demonstrated. 

PACS numbers: 03.65.Bz; 03.67.-a; 03.67.Hk 



I. INTRODUCTION 



The discrimination of quantum states is a nontrivial 
problem since a quantum state cannot be cloned per- 
fectly if it is unknown to us There are strategies of 
reaching the optimal results in the discrimination mea- 
surements: the minimum error probability discrimina- 
tion 0, HJ , where each measurement outcome select one 
of the possible states and the error probability is mini- 
mized, and the optimal unambiguous discrimination for 
the linearly independent states [3, H, IE 0> HI > where un- 
ambiguity is paid by the possibility of getting inconclu- 
sive results from the measurement. Recently an analogue 
of the optimal unambiguous discrimination is proposed 
for the linearly dependent state sets as the maximum 
confidence measurement Q . In all these approaches, the 
input states to be discriminated occur or are prepared 
with some a prior probabilities in reality. 

If the states we want to discriminate are known, we 
can use the given information to find the optimal mea- 
surements, which can be either von Neumann measure or 
the generalized quantum measure in the form of positive 
operator value measure (POVM), to discriminate the ele- 
ments in the set of states. If we know nothing about these 
states, however, the only information available will be the 
permutation symmetry with respect to a given number 
of copies of the unknown states provided to us. In the 
original work of two completely unknown states discrimi- 
nation without ambiguity [lfj , the authors used a sort of 
programmable quantum devices, which have been stud- 
ied both theoretically and experimentally in the recent 
years [H QJ, 0, M, QJ, [H, [13, O , to relate the pro- 
gram part of the inputs in a simple way to the unknown 
qubits IV'i) and \ip2) that one is trying to identify. The to- 
tal input states measured by the device are thus prepared 



with the pairs of unknown qubits as |?/>i) aIV'i) b\^2) c an d 
IV'i) a\^2) b\^2) C i a nd the optimal POVM for the unam- 
biguous discrimination of the inputs is designed with the 
permutation symmetry of the program registers A, C and 
the data register B. If the a priori probabilities of the in- 
puts are equal, the maximum average success probability 
of discriminating them can be as large as 1/6. 

The most general problem of this type is when we have 
riA copies of the state in the program system A, nc copies 
of the state in the other program system C, and tib copies 
of the state in the data system B. Then the task is to 
discriminate two input states 
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(1.1) 
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with the minimum error or with the least inconclusive 
probabilities, if we apply the minimum-error or the un- 
ambiguous strategies, respectively. We should optimally 
distinguish between the above inputs with respect to 
these cost functions, keeping in mind that one has no 
knowledge of and \ip2}- The optimal schemes, in 
which the multiple copies of program or data used in the 
inputs, were solved for the input states prepared with 
equal a priori probabilities but with an arbitrary num- 
ber of copies in the program registers (ua = nc = n, 
iiB = 1) 0, for the inputs with multiple copies 
of data (riA — nc = 1, Ub — w)[21[, and for the 
unambiguous state discriminator working in the whole 
range of the a priori probabilities of the inputs with 
jia = nc = n and tib = 1[22||. The programmable 
discriminator to unambiguously discriminate a pair of 
unknown input states prepared with single program and 
data copies (jia = ns = nc = 1) is also generalized in 
the minimax approach [23| to that of N states [24{ . 

In this paper we study the unknown qubits discrimi- 
nation with the inputs prepared with n program copies 
and m data copies (ua = nc = n, jib = rn). Since 
the efficiency in discriminating the averages of the inputs 
in Eq. (1.1) indicates the upper bounds we aim to ap- 
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proach in the discrimination of a pair of unknown states, 
we will study the problem of discriminating the averaged 
input states as in [H, [2(| [2l[ by both unambiguous and 
minimum-error discrimination strategies and generalize 
part of the results in these papers. We will also show the 
dependence of the optimal solutions on the a prior proba- 
bilities of the averaged input states. To apply the Jordan 
basis method [25| to the optimal discriminations of the 
mean input states < |$j n )(*i n | > and < >, 
we give a systematic study of the structure of these mean 
input states through finding an inherent symmetry which 
exists only under the condition ua — nc- In the limit 
of very large numbers of both program and data copies, 
we demonstrate that the definite discrimination with a 
unit success probability for the average input states will 
be approached for all a priori probabilities. 

This paper is organized as follows: in Section II we 
demonstrate the equivalence between the discrimination 
of the average input states and that of uniformly dis- 
tributed mixed states. The structure of the average input 

I 



where C% is the number of ways to choose k objects from 
a group of n objects without regard to order. The index k 
of |efe) means how many digits 1 in this basis vector out of 
the total digits, and these basis vectors are orthonormal, 

{&i | &j ) ' 

As in the related works, we assume that the unknown 
qubits IV'i) and ^2) uniformly and independently dis- 
tribute over the Bloch spheres, and the averages of the 
inputs with ua — nc — n and ub = m are given as 

91 = j4^J d ^J d ^" l )(n n \ 



mixed states is discussed in Section III and Appendix A 
and B, and the derivations of the inner products of the 
Jordan basis, which are used to represent the mean in- 
put states, and their multiplicities are given in Section 
IV. Our main results about the optimal unambiguous 
discrimination and the minimum-error discrimination of 
the average input states are given in Section V and VI, 
respectively. Finally, we give some conclusive discussion 
in the last section. 



II. MEAN INPUTS AS UNIFORMLY 
DISTRIBUTED MIXED STATES 

In this section we demonstrate that the discrimination 
of the average input states is equivalent to that of two 
uniformly distributed mixed states. First we use the bi- 
nomial expansion of n copy tensor product of qubits, 



(2.1) 



(2.2) 



(n+l)(n+m+l) 

= („ + l)(n + m + l) E MM 

(n+l)(n+m+l) 

- (n + l)(n + ro+ l) £ MM. (2-3) 

v i— 1 



where \vi) = \ej) a.b ® \ e k)c {0<j<n + m,0<k<n) 
and \v[) = \ej) A ® \ek)B,c (0 < j < n, < k < n + m), 



(cos(0/2)|O> +sin(0/2)e^|l))® n = ^ cos"- fc (0/2) sin fe (0/2)e 1 ^ J C*\e k ) , 



fc=0 



to introduce in the orthonormal basis {|ej)} 

1 



\ei) 
\e k ) 



1 
1 



|0,0,...,0) 

(j 1, 0, . . . , o> + jo, 1, - • - , o> H — + jo, 0, . . . , 1>) 

( 11,1,..., 0,0) +|0, 1, 1, . . . , , 0) + • ■ ■ + |0, 0, . . . , 1, 1» 

k' s 1 in n digits 

S v ' 

summation of terms 



'Cm 



|1,1,...,1) , 



3 



and we have also used the integral, 

2 r sin 2 ™-^ cm*- 1 ! dx = r(m)r(n j . (2.4) 
Jo r ( m + n ) 

In taking these averages we actually realize the one- 
to-one maps from the unknown qubit ensembles {|^>i)}, 
{l^a)} (the wave bracket means a set) to two known 
mixed states: 

— » pi 

{|V>2>} — > Pi. (2.5) 

Generally these mixed states are produced from the un- 
known qubit ensembles with some different a prior prob- 
abilities rj and 1 — rj, respectively. In the following we 
will apply the Jordan basis method [25] to derive the 
optimal schemes for the unambiguous and the minimum- 
error discrimination of the mixed states p\ and p2 ■ 

III. STRUCTURE OF MEAN INPUT STATES 
A. Closed Chains of Basis Vectors 

We here study the structure of the mean input states 
pi and p2- Let H\ be the Hubert space of p±, which 



is spanned by and H2 the Hilbert space of P2, 

which is spanned by {W)}. We have dimHi n H 2 = 
2n + m + 1 , and the dimension of the total Hilbert space 
H is therefore 



dimH = dimHi U H 2 = dimHi + c?im7J 2 — dimHi D H 2 
= 2n 2 + 2nm + 2n + m + l. (3.1) 



For a particular \vi) in Hi, we need to find out with 
which elements of {|t>l-)} spanning H2 it has non-zero 
overlaps , i.e. to find out all \v'j)'a such that (v$|t^) ^ 0. 
To do it, we use the fact, 



where i+j = k, to split \e.j)A,B in \vi) and \ek)B,c i n 
\v[) into the summation of tensor products by parts as 
follows: 



|efe)A,i 



\ei)A\ej) B + 



jf-fi+l/^j—l 
i On O m 



|e»+iM|ej-i)s 



(3.3) 



where i + j = k, and the number I is determined by how 
many ways two non-negative integers can be summed up 
to be k. 

Next, we use the [714, ri2, W3] symbol defined as follows 



to classify the mutually overlapped subspaces of Hi and 
H2- With the help of the above formula, we find that 
any basis vector \v p ) in Hi or \v' p ) in H2 has a unique 
representation in terms of this symbol: 



(Ji (J3 

\h^r\ e ^ A \ e o)B\^k)c + 

V ^n+m 

[i,j,k] + [i + lj-l,k} + - 



/ On O m 

+ [i + l,j-l,k], 



e i+ i)A|e ? -_i)s|e fe )c 



i+l/tj — l 



' On O f 



-| e i+/>A|ej-/)s|e fe )c 



(3.4) 




e/b)A|ei)s|ej)c 



' On Om 



efc)A|ei+i)s|e J -_i) c 



' On Om 
^n+m 



|efc)A|e»+i)s|ei-j)c 



(3.5) 



where the coefficients of in each terms of the basis vec- 



tors have been absorbed into the square brackets. Obvi- 
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ously we see from this expression that any couple of \v p ) 
and | Up) satisfy (v' p \v p } + 1 only if they have a common 
[711,712,713] term. 

By the [711,712,713] representation of \v p ) or \v' p ), each 
term of basis vector has the same number, N = n\ + 112 + 
713, which goes from to 2n + to. We'll show that each 
distinct N corresponds to a pair of closed basis vector 
chains, one of which is in Hi and the other of which in 
if 2, and their elements may have non-zero overlaps. 

Let's first look at N = and N = 2n + to, the two 
simplest cases. [0,0,0] = |0, 0, • • • , 0) and [n,m,n] — 
|1, 1, • • • , 1) are shared by Hi and H2, and their overlaps 
are just (v'\v) — 1. They are two initial pairs of closed 
basis vector chains with only one element. 

When N = 1, let's pick out one basis, [0, 1, 0] + [1, 0, 0], 
in Hi (by [711,712,713] representation all the terms of a 
basis in Hi have the same last digit 713), and then one 
basis, [0, 1, 0] + [0, 0, 1], in H2 (the first digit »i are the 
same in [ni,n 2 ,n 3 ] representation), which is overlapped 
with it. The only other basis vector in Hi, with which 
[0, 1, 0] + [0, 0, 1] of H 2 also has overlap, is [0, 0, 1], while 
[0, 1, 0] + [1, 0, 0] in Hi is also overlapped with [1, 0, 0] in 
Hi- Thus we exhausted all basis vectors in Hi and H2 
with N = 1 and obtain such a pair of closed chains as 
follows: 

|«i) = [0,1,0] + [1,0,0], 

\v 2 ) = [0,0,1]; (3.6) 

K) = [0,1,0] + [0,0,1], 

14) = [1,0,0]. (3.7) 

Both of them have 2 elements and, by retriving the co- 
efficients absorbed in the square brackets, we can easily 
find all their overlaps: 



(4h> = Mv£> = a/tP-. (3-8) 

For N — 2n + m — 1, we find in the same way the 
following two closed chains: 

\vi) — [n, m — 1, n] + [n — 1, to, n], 

\v 2 ) = [n,m,n- 1]; (3.9) 

\v'i) = [n, m — 1, n] + [n, m, n — 1], 

|4) = [n-l,m,n]; (3.10) 

and, with the combinatorics identities, we find that they 
have the same overlaps as in N — 1 case, so we call them 
the conjugate chains of N = 1. 

As N increases to any i < n, the closed chains of basis 
vectors can be found inductively in the above way: 

= [0, t, 0] + [1, <- 1, 0] + [2, i- 2,0] + ••• + [*, 0,0], 



\v 2 ) = [0,7 -1,1] + [1,7 -2,1] + ••■ + [1-1,0,1], 

|«s) = [0,t- 2,2] +••• + [* -2,0,2] 
\v i+ i) = [0,0,i]; (3.11) 

\v[) = [0,t,0] + [0,*-l,l] + [0,*-2,2] +.-. + [0,0,*], 
|4) = [1, 7-1, 0] + [l,7-2,l] + --. + [l, 0,1-1], 
|4) = [2,t-2,0]+.» + [2,0,7-2] 

K+i) - M,0]. (3.12) 

For every \vj) (J < i+ 1) in a chain of Hi, the last index 
713 = j — 1 increases as j from to i, while for every \v'j) 
(J 5: i + 1) in a chain of H2, the first index n 1 = j — 1 
increases as j from to i. The number of different terms 
in a (resp. \v'j)) is how many ways the non-negative 
integers 711 and 712 (resp. 712 and 72,3) can be added up to 
i — j + 1, so it goes down from i + 1 to 1 as j increases. 

If iV = 2n + ?7i — i, there are two chains conjugate to 
those of N = i: the last index 713 of \vj) decreases from 
n to n — i as j goes up, while the first index 711 of \v'j) 
also decreases from n to n — i as j goes up. The number 
of different terms in a \vj) (resp. |4)) is how many ways 
the non-negative integers ni and 712 (resp. 712 and 713) 
can be added up to n + m — i + j — 1. 

In the pairs of closed chains with mutual overlaps, 
whenever N increases its values by 1, the number of 
the elements in the closed chains will increase by 1, if 
< N < n. On the other hand, the number of the ele- 
ments in closed chains also increases by 1 as N decreases 
by 1, if n + m < JV < 2n + m. However, the last index 713 
of the \vi) in Hi and the first index ni of the |t^) in H 2 
cannot increase beyond n or decrease beyond 0, so the 
number of elements in a closed chain cannot be larger 
than 77 + 1. 

The size of a closed chain thus increases in two direc- 
tions from both N — and N — 2n + to, and it will be 
fixed with n + 1 elements as N increases to the medium 
value 77 or decreases to the medium value n + m. The 
number of such chains with the maximum n + 1 elements 
in Hi and H2 is to + 1. Therefore, the total number of 
the elements in all closed chains in Hi or H2 equals 

n 

dimHi = dimH2 — 2 i + (to + l)(n + 1) 
i=i 

= (71 + TO+ 1)(t7 + 1). (3.13) 

After we have obtained all pairs of closed chains with 
mutual overlaps this way, we will find that the {|ui)} in 
a closed chain of Hi and its counterpart {|4)} m -^2 can 
be permutated such that they have symmetric overlap, 
( v i\ v j) = ( v i\ v j)> f° r an y couple of i and j (we can also 
have i = j because the inner products are real). We call 
it a mirror symmetry , and give a procedure of how to 
obtain pairs of closed chains with this symmetry in Ap- 
pendix A. This symmetry leads to very useful invariants 
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for finding the Jordan basis inner products for each pair 
of closed basis vector chains. 



B. Mean Input States Represented by the Jordan 
Basis 



As it is proved for any pair of subspaces j26[, there 
exist the Jordan basis {|</>i)} in Hi and {]</>£)} in H2 , 
with which the mean input states p\ and p 2 are expressed 
as 



Pi 



f>2 



(n + l)( 



(n+l)(ra+m+l) 

l)(n + to + 1) ^ 

(n+l)(n+m+l) 

TT7 TT y \<t>i){<f>il (3-14) 

(n + l)(n + m + l) ^ 



i=l 



where (<^j|^-) = for each couple of i ^ j. If we have 
obtained the mean states with the above representation, 
we will be able to use the strategy in Ref. [25[ to get 
the optimal scheme for the unambiguous discrimination 
of them. 

To get the mean inputs represented by the Jordan ba- 
sis, we need to look for the orthogonal transformations, 



/ I'M \ 




( an 


a-12 ■ 




1 


\vx)\ 








a-21 


a-22 ■ 


• &2k 




\V2) 


(3.15) 


V \M J 




\ Clkl 


ak2 ■ 


■ O-kk J 


V 


\Vk) j 





in Hi, and the corresponding orthogonal transforma- 
tions, 



/ Wi) \ 

W2) 

\ Wk) J 



( 



\ 



Hi 
'21 



'12 

'22 



'1.2 



'Ik 

l 2k 





( M> \ 




W2) 


) 


K K) ) 



(3.16) 



in H2, for all pairs of closed basis vector chains (the in- 
dex k means the size of a closed chain). Though it is 
very difficult to obtain all these orthogonal transforma- 
tions analytically, in the following we will use the sym- 
metric property of the basis vector chains to derive all 
the Jordan basis inner products, for each closed 

basis vector chain, which are enough for us to find the 
optimal POVMs and their success probabilties for the 
unambiguous and the minimum-error discriminations of 
our average input states. 

By the inverse orthogonal transformations, \vi) = 

EfcOwl^O and K) = E fc a*M>. of Eqs- (3.15)-(3.16), 
we deduce from the mirror symmetry, (v'^Vj) = (vi\v'j), 
the following equations, 



k 



a ik a 



E 

k 



a'ik^j^'kl^k), (3-17) 



for each pair of i = j or i ^ j. We have the sufficient con- 
dition, a,ij — a'ij for all i's and j's, to guarantee the valid- 
ity of all these equations and, therefore, two orthogonal 
transformations in Eqs. (3.15)-(3.16) can be reduced to 
a single one in realizing the transformation to the Jordan 
basis. 

Substituting |^) = J2k a ik\ v k) and |$) = Y,k a ik\ v 'k) 
into Ei^l^i) and using the orthogonal transformation 
property Ei a fc'i a «fc = ^fefe'j we wm fi n( l the following 
invariants, 

N N / \ 

s N = e<«> = E (Ev^k'i)(E a ^» 

i=0 i=0 V k' k ) 

N 

= 5>iK>, (3.18) 



i=0 

if < N < n— 1. For the conjugate chains with n + m + 
1 < N' < 2n+m we have Sn> = S2n+m-N, and the m+1 
pairs of closed chains with the maximum n + 1 elements 
all have the same invariant. These invariants are very 
useful in finding all the Jordan basis inner products, and 
they are determined by how many elements in a closed 
basis vector chain. We give a simple method of how to 
calculate these invariants with the [ni, 712,713] symbol in 
Appendix B. 



IV. DERIVATION OF JORDAN BASIS INNER 
PRODUCTS AND THEIR MULTIPLICITIES 

For N — - and N = 2n + to the Jordan basis are 
just the original ones, and the inner product is obviously 
(4>i\4>[) = 1, i.e. they belong to the intersection of two 
subspaces. 

In N = 1 (similarly in N = 2n + m — 1) case, we 
apply the rotation Ti 2 on the closed chains given by Eqs 
(3.6)-(3.7), 



>2> 



= Zl, 2 

= T h2 



l«2> 



cos 9 sin 9 
— sin 9 cos 9 



\V2) 



cos tt 
— sin( 



sin t 
cos ( 



•(4.1) 



to obtain l^a)} and {\(/>'i), \ 4Q}- The conditions 

(^|0 2 ) = (0i|0 2 )=O, (4.2) 

imply 

- (vi\v[) tan 6» + (vi\v 2 )(l - tan 2 9) 



0, 



(4.3) 



a quadratic equation that can be solved easily. Substi- 
tuting tan 9 and (filt^) m t° the inner products 
of Jordan basis, we obtain 



1. 



(4.4) 
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We can verify that their summation equals to the invari- 
ant Si, which is the square of the coefficient of [0,1,0] 

c 1 

term, jts^- 

The TV — 2 Jordan basis inner products are obtained 
by first applying two successive rotations T = Ti^Ti^ in 
3 dimensional subspace and then imposing the vanishing 
inner products for the pairs of rotated vectors with dif- 
ferent indices to seperate one couple of basis vectors from 
others; after that we use T 2i 3 to seperate the rest pairs 
in the same way. The rotation Tij in the subspace of a 
chain is defined asannxn(n is the number of chain ele- 
ments) identity matrix with the elements in the positions 
and {j,j} replaced by the correspond- 
ing elements of an 0(2) matrix. Following the procedure, 
we obtain the inner products with their absolute values 
arranged in a descending order: 



(0i|0i) 

(02 |0' 2 ) 



n + m 

n(n — 1) 
(n + m)(n + m — 1) ' 



(4.5) 



The summation of these Jordan basis inner products is 
equal to the invariant S 2 , which can be independently 
obtained by summing up the square of the [0, 2, 0] coeffi- 

cient, ci m , and the square of [1,0, 1] coefficient, c i " , 

if m > 2. If m < 2 instead, we use only [1, 0, 1] term to 
find the same result for m = 1. 

We see from the above results that the first two inner 
products of TV — 2 chains are those of TV = 1 chains, and 
the first inner product is just that of TV = chains, so 
this recurrence gives rise to the multiplicities of Jordan 
basis inner products in the whole Hilbcrt space. 

For any pair of closed chain with a fixed TV num- 
ber, if we apply successive orthogonal transformations, 

J-l — J-1,NJ-1,N-1 ■■■ J- 1,3-11,2, ^2 = -£2,AT-£2,AT-1 • • • -t2,3, 

• • •, and impose the vanishing inner products for the pairs 
of basis vectors with different indices, we will obtain the 
Jordan basis inner products as in the above equations. 
A useful obsevation in doing so is the recurrence of the 
Jordan basis inner products from closed chains to closed 
chains, which is found by induction as we have done for 
the chains of TV = 0, 1, 2. The general recurrent pattern 
for these Jordan basis inner products is as follows: the 
inner products in the chains of TV = k coincide with those 
of TV = k + 1 chains except for the one extra inner prod- 
uct of TV = k + 1, if < k < n — 1; meanwhile, the 
inner products of TV — k chains are those of TV = k — 1 
except for the one extra inner product of TV = k — 1, if 
n + m+1 < k < 2n + m. 

Together with the invariants Sn, we can tremendously 
simplify the calculation to get the Jordan basis inner 
products with this recurrence. If we have obtained the 
Jordan basis inner products for the chains of TV = k 
(0 < k < n — 1), the one more unsolved inner product 
for the chains of TV = k + 1 is just Sk+i — Sk, which can 



be independently obtained from the closed chains repre- 
sented by the original basis {|t>i)} and {|t^)}. Starting 
from TV = 0, we can find all the Jordan basis inner prod- 
ucts for the chains of0<TV<n — 1 this way. For 
example, the extra Jordan basis inner product for the 
TV = 3 chains is obtained by S3 and Eq. (4.5) as follows: 



(04|44) 



= s 3 

Cm 



S 2 



+ 



S~1l s~tl 
f<2 

n(n 



cl 



cl 



C 2 
l)(n-2) 



(n + m)(n + m — l)(n + m — 2) 



, (4.6) 



where S3 is from the contribution of [0, 3, 0] and [1,1,1] 
terms, if we suppose m > 3 (otherwise we can check out 
the above result too with the relevant terms for m = 1 
and m = 2, respectively). For the chains oi n + m + 1 < 
TV < 2n + m, we can use the similar procedure to get 
all their Jordan basis inner products too. Therefore, the 
m + 1 pairs of closed chains with the constant invariants 
Sat for n < TV < n + m have the following complete set 
of n + 1 Jordan basis products: 

(01 101) = 1, 

n 

(0 2 |0 2 ) = 



(03|0 3 ) 

(0™+iK+i) 



n + m 

n(n — 1) 
(n + m)(n + m — 1) ' 

n(n — 1) • • • 1 



± 



(n + m)(n + m — 1) • • • (m + 1) 



TV «"> 



The sign for the last inner product is determined by 
whether n is even or odd. Thus we have exhausted all 
the closed basis chains and obtained all their Jordan basis 
products. 

As the conclusion for this section, we list all inner prod- 
ucts of Jordan basis and their multiplicities in the follow- 
ing table: 



inner product 


multiplicity 


1 


2n+m+l 


n 
n+m 


2n + m — 1 


n{n— 1) 
(n+m) (n+m— 1) 


2n + m - 3 


n(7i-l)(n-2) 


2n + m — 5 


(n+m)(n+m— l)(n+m— 2) 






n{n— 


m + 1 


(n+m)(n+m— l)---(m+l) 



The inner products, (0|0') = 1, correspond to the inter- 
section of Hi and H 2 , and we also see that each closed 
chain in Hi has a one-dimensional joint space with the 
corresponding chain in H 2 - 
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V. OPTIMAL UNAMBIGUOUS 
DISCRIMINATION OF AVERAGE INPUT 
STATES 

In Ref. [25|, the optimal scheme for the unambiguous 
discrimination of two mixed states represented by the 
Jordan basis, 



i 

P2 = I 



(5.1) 



which occur with the a prior probability of rj and 1 — rj, 
respectively, is derived through optimizing the POVM 
in the form of IIi = ^iLi n i, 4 , U 2 = J2i=i n 2,» an d 
n = / — III — n 2 with / being the identity operator, 
where 



and 



IIi,, 



l - q% 



1-K^>l a 
1-mm 2 



\Zi){Zi\, 



\Vi)(Vi 



(5.2) 



(5.3) 



The orthonormal basis and in the above 

equations satisfy 



l(«l 2 = i-m\<k)\ 3 , 



i\Vi)=0, (5.4) 



for all i's and j's, and qt, q\ the failure probabilities in 
the unambiguous discrimination of couples of \<fii) and 
|$) for all i's. Then the total failure probability for the 
discrimination of two mixed states, 



i=l 



is optimized with Qi taking the following solutions: 



yotj + (1 ~ 7 l)Pi\ { ( t ) 'i\ < t ) i)\ 2 ifv<Ci 

QT = { 2y/r}(l-v)<XiPiM\<f>i)\ if c * < V < * , 



r,aim\(t>i)\ 2 + (l-v)f3i if V>di 



(5.6) 



where Cj, di are the boundaries of the range given as 



L = \ci,di\ = 



(5.7) 

In the range of Jo = C\ih, where the total POVM is valid, 
the optimal POVM taking all absolute minimums on the 
center line of Eq. (5.6) achieves the success probability, 



p(n) = 1 - qt = 1 - 2^iTT)J2 V^mmi (5- 

i 

if the spaces of two mixed sates to be discriminated join 
at the origin. 

In our problem, however, there is a non-trivial inter- 
section space of the Hilbert spaces H\ and H2, so this 
intersection contributes to failure probability and should 
be added to the total failure probability Ql too. With 
all the Jordan basis inner products and their multiplic- 
ities derived in the last section, we can readily obtain 
the optimal scheme for the unambiguous discrimination 
of our average input states prepared with any number of 
program and data. 

Taking the on and /3j in Eq. (3.14) in the general sit- 
uation with the input states produced with n copies of 
program and to copies of data, we find the validity range, 



(n + to) 



(n + to) 2 + n 2 ' [n + m) 2 + n 2 



(5.9) 



^liVQHQiiv) + (1 — r /)Piq'i( r l)} ; (5-5) and the optimal success probabilty, 



op l 



2n + to + 1 1— 1 n(2n + m — 1) n(n — l)(2n + m — 3) 

1 — — — 2 \/ 77(1 — rj)-. — r{ 1 — — — 

(77, + to + l)(n + 1) [n + m + l)(n + 1) n + m [n + m)(n + m — 1) 

77,(77, — 1)(tt, — 2)(2n + to — 5) n(n — 1) • ■ ■ 1 



(n + m)(n + to — l)(n + to — 2) [n + m){n + to — 1) • • • (to + 2) 

of the total POVM. We have 

2n + to + 1 2771 + 2 



}, (5.10) 



P(rj) < 1 - — — - — - — x if(n,m)vw - (5-H) 

w/ (n + TO+ l)(n+ 1) (n + to + l)(n+ 1) v JV x 
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and 



where 



2n + to + 1 An + 2m - 2 



> 1 - f T in7 T^n - ? I mv ZTi x tffa^vW^, (5-12) 
(n + m + l)(n + lj (n + to + l)(n + 1) 



n n(n — 1) n(n— l)(n — 2) n(n— 
K(n, m = + y -, '- + 1 + ■ • • + -, ^ '— -. 5.13 

n + m (n + m)(n + to — 1) [n + m){n + m — ±){n + m — 2) (n + m){n + m — 1J • • • (to + 1) 



is a typical convergent series. From these inequalities, we 
see that P{rj) will tend to a certain 1 if both n and m 
tend to infinity. 

We here give some examples when the optimal success 
probabilities P{rj) can be reduced to the closed forms: 

For one copy program and to copies data in the input 
states, we have 



is solved by classifing the enginevalue spectrum {Xk} of 
the operator Q, 



dimH 



A = 1]2P2 - VlPl = J! ^k\u k )(uk\- 



(6.2) 



fc=i 



P(V) = i - 



1 



2(m + 2) to + 2 



x y/vO- ~ vY (5-14) 



For 2 copies of program and to copies of data used in the 
input states, we get 



P(V) 



1 



to + 5 
3(to + 3) 

4(to + 4) i— , 

3(to + 2)(to + 3) X (5 ' 15) 



For n copies program and 1 copy data in the input states, 
the optimal success probability is found as follows: 



p( v ) = l- 



n 
1 



x 207(1 - V) 



x- 



n 
1 - 



2 (n+l)(»+2) 
(n 2 + (n - l) 2 + • • • + 2 2 + l 2 ) 



The corresponding optimal measurement is achieved by 
two projectors, 

feo-l 

rii = \uj k }(uj k \ 

k=l 
dimH 

n 2 = Yl kfc)<Wfc|, (6.3) 

where III is the projector onto the space spanned by 
fco — 1 eigenvectors of A with negative eigenvalues and 
II2 the projector onto the space spanned by the remain- 
ing eigenvectors with non-negative eigenvalues, and they 
satisfy IIi + II2 = /, where / is the identity operator. 
The minimum error probability is therefore given as 

Pe = ~(1 - Tr|A|) = i(l - Tr\ mP2 - mPl \), (6.4) 



2n(2n+l) r — . 

3(n + l)(n + 2) 



n + 2 

If 7/ = 0.5 (equal a prior probabilities), it reduces to 

n 



with |A| = x/AU. 

With the Jordan basis representation of our mean in- 
put states, the corresponding A operator is given as 



A = 



1 



(n+l)(n+m+l) 



P 



3(n + l) 



(5.17) 



(n + l)(n + to + 1) 



i=l 



given in [2(| and tends to 1/3, the optimal IDP average 
[J, [H, [f| in discriminating a pair of known qubits, as n 
goes to infinity. 



VI. MINIMUM-ERROR DISCRIMINATION OF 
AVERAGE INPUT STATES 

The general problem of discriminating two mixed 
states pi and pi , which occur with different a prior prob- 
abilities T]i and 7/2 , respectively, with the minimum of the 
error probability, 



where Aj = 772 1 0^) ($| — vA't'i) {4*^- We derive its eigen- 
value spectrum in the following. 

The subspace spanned by each pair of {|&),|<^)}, 
where 1 < i < (71 + l)(n + to + 1), is orthogonal to 
the rest part of Hilbert space H. In it we introduce in 
the following new bases: 



K-) = \4>i) 



(6.6) 



if = 1; and 



1 



Pe = 7?iTr(pin 2 ) + rfrTripzUi), 



(6.1) 



v/2(l + <&|#» 
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if (4>i\4>i) + 1 1. All these new bases are orthonormal 
((oJi\uJj) = S hj , (cj'^j) = Sij and (Wj|u^) = 0). 
With this new basis, the Aj in the second case is given 

as 



1 { miam) M\a+M) 



\(m - t? i)(i + Kj) 



, (6-8) 



where k, = 



We thus obtain the following eigen- 



values of the A,- 



W = ^(c+Vl-(l-c>?) 



(6.9) 



where c = 772 — 771 . The eigenvalue spectrum of A is 
therefore obtained as follows: 



2n+m+l (n+m+l)(n+l) , . 

(n + l)(n + m + l) ^ 2(n + l)(n + m + l) l=2r f^ n+2 V v / 
( n+m+ i)(„+i) 

2(n + l)(n + m + l) i=2 ^ m+2 V V J 

I 



where and |AJ) are the eigen- vectors corresponding 
to the eigenvalues and respectively. 

Together with the table of Jordan basis inner products 
((f>i\4>i) in Section IV, this eigenvalue spectrum of A allows 
us to obtain the minimum error probability for arbitrary 
n and m in the input states. Here we give two examples 
when c = (equal preparation probabilities 7/1 = 7/2) and 
compare their results with those of optimal unambiguous 
discrimination we obtained previously. 

First, we take 1 copy of program and n copies of data. 
In this case there is a multiplicity of n + 1 for (4>i\4>' i ) = 
— l/(n+ 1) and none of other Jordan basis inner product 
is not equal to 1. Plugging these results into (6.4), we 
obtain 

As 71 goes to infinity, Pe tends to 1/4. In the unambigu- 
ous discrimination for this case, the least failure proba- 
bility Ql has the limit of 1/2 as n tends to infinity. So 
we have the relation, Pe = 0.5Ql, if we have 1 copy of 
program and infinite copies of data in our inputs. 

Another example is to have 71 copies of program and 
1 copy data. After substituting the Jordan basis in- 
ner products and the corresponding multiplicites into the 
Helstrom formula Eq. (6.4), we obtain the minimum er- 
ror probability 

P £ = l(i-_-_yVi-(_--)2 — V (6.12) 

E 2y n + 2^V y n+l' n+lj 1 ' 



If n — ► 00, the limit of Pe is 

P E = X - (\ - 2^ y/l-x 2 x dx^j = i (6.13) 

which is consistent with the result obtained by other 
method in [l9( together with those of all other 77. Com- 
pared with the corresponding least failure probability, 
Ql = 2/3, in the unambiguous discrimination, we have 
Pe = \Ql if we are dealing with the inputs carrying 
infinite copies of program and 1 copy of data. 

The general relation [28| . Pe < \Ql, for the error 
probability Pe and the failure probability Ql of the 
unambiguous state discrimination is well satisfied for 
these extreme situations. Moreover, we see from the re- 
sults a prominent difference of the minimum-error dis- 
crimination from the unambiguous discrimination: more 
copies of program copies will give higher success probabil- 
ity, while in the unambiguous discrimination more data 
copies gives higher success probability. 



VII. CONCLUDING REMARKS 

So far we have studied the problem of discriminat- 
ing pairs of unknow qubits with the devices processing 
the inputs prepared with n and m copies of them used 
for the programs and the data, respectively. The upper 
bound efficiency of the devices are achieved by optimally 
discriminating two uniformly distributed mixed states if 
we assume the pair of qubits we want to discriminate 
uniformly and independently distribute over their Bloch 
spheres. We show that there exists a symmetry in these 
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mixed states when number of state copies in two pro- 
gram registers are equal. If we represent these mixed 
states with the Jordan basis, we can conveniently obtain 
the Jordan basis inner products and their correspond- 
ing multiplicities with this symmetry. Then we have the 
optimal schemes of the unambiguous and the minimum- 
error discrimination of these mean inputs occuring with 
arbitrary a prior probabilities. 

We should also compare the upper bound efficiency of 
the unknown state discriminators we discussed with that 
in discriminating a pair of known pure states |^>i) and 
|^2) • Given m copies of such known states, the upper 
bound success probability |, S 1, 1 - \(i/>i\ifa)\ m , for 
their unambiguous discrimination, which is lower than 
that of the minimum-error discrimination, promises a 
definite identification if m tends to infinity. Correspond- 
ing to our unknown input states with m copies of data 
in Eq. (1.1), however, we see that the definite discrimi- 
nation can be approached only if both n and m are very 
large numbers. In the design of any programmable dis- 
criminator for unknown states, we will need to use large 
program registers to improve its efficiency. 

For the situation of arbitrary tia , ns and nc, the mir- 
ror symmetry we find in our mean input states is gone, so 
it is impossible to obtain the analytical solutions to the 
optimal discrimination problems as what are discussed in 
this paper. In this even more general situation we need to 
deal with the mixed states of different dimensions. With 
the structure of the mean input states we clarify in this 
paper, however, a numerical algorithm can be worked out 
to obtain their Jordan basis inner products, since there 
exist the Jordan basis for a pair of arbitrary subspaces 
[2c| . It is also possible to generalize the optimal solu- 
tions to those of unknown qudits discrimination. These 
problems will be studied in the future work. 



APPENDIX A: PROCEDURE TO OBTAIN 
SYMMETRIC CLOSED BASIS VECTOR CHAINS 

Here we give a procedure of obtaining the closed ba- 
sis vector chains with the mirror symmetry, (Uj|i>»-) = 
(v'ilvj), for any couple of indices i and j. 

Step 1: Generating the closed basis vector chains of 
any fixed N = m + n 2 + n 3 . 

In Hi , since each term of a basis vector in this repre- 
sentation has the same digit n 3 , we start from n 3 = and 
obtain the corresponding basis vector (that is to find out 
all the combinations of non-zero intergers that added up 
to be N). Then we increase n 3 by 1 each time till 713 = n 
and obtain in this way all basis vectors in this closed 
chain. Similarly we obtain the corresponding chain in 
H 2 , each element of which has the same first digit n\ for 
all terms. 

Step 2: Permutation of the basis vectors, which is given 
by [m, ri2, TI3] representation, in the obtained closed basis 
chains. 



Remember that k] 



1 \ 



\ei)A\ej)B\ek)c m 



Hi and [k,j,i] 



\ek)A\e-j)B\ei)c hi H 2 . Pick 



out in a closed chain all couples of basis vectors carrying 
[Zi , Z2 , ^3] and [Z3, Z2 , /1] terms (the terms with n\ and n 3 
interchanged), which are guranteed to exist by the fact 
that N is constant for any of a closed chain, and let the 
basis vector with [1%, l 2 , l 3 ] be \vt), and that with [Z3, l 2 , h] 
be |t!j_|_fc), where i, k are the integers within the range of 
the size for a closed chain, in H\. Correspondingly in H 2 , 
we set the basis vector carrying [l 3 , l 2 , h] instead be 
and that carrying [h,l 2 ,h] be \v' l+k ). 

Then we proceed from other terms in the basis vectors 
to obtain two whole chains with the mirror symmetry. If 
the number of program copies ua and nc are unequal, 
we cannot guarantee the existence of the couples in the 
form {[Zi, l 2 , 13], [l 3 , l 2 i £1]}) and the mirror symmetry of 
the basis vectors doesn't exist. 

With a not so large TV, we here give an example of how 
to perform the procedure. For the case of 4 copies of 
program and 1 copy of data in the input states, suppose 
that we need to obtain the chains of N = 4 with the 
mirror symmetry. 

There are 5 elements in the chains as we studied previ- 
ously. Starting from n 3 = and n\ = 0, we see that 
the couple {[h, l 2 , h], [h, h, h]} in the basis vectors is 
{[4,0,0], [0,0,4]}. Adding the rest terms in these basis 
vectors, we set \vi) = [4, 0, 0] + [3, 1, 0], \v 2 ) = [0, 0, 4] in 



Hi, and K) = [0,0,4] + [0,1,3], 



[4,0,0] in H 2 . 



These two couples of basis satisfy (vi\v 2 } = (i>2|i>i) = 



c°cl 



The basis in H 2 carrying [3,1,0] is [3,0,1] + [3,1,0], 
and we set it as \v' 3 ). In Hi the basis carrying [0, 1,3] is 
[1,0,3] + [0, 1,3], which is set as \v 3 ). We have {vi\v' 3 } = 

Then we set |u 4 ) = [3,0, 1] + 



\V-l\V3) -- 

[2,1,1] and 

(^ 3 M = 



~cT 



cic\ 



v.. 



= [1, 0, 3] + [1, 1,2], and there is (v 3 \v' 4 ) 



C}C° 



ci 



clcl 



Finally we consider the terms with [2, 1, 1] and [1, 1, 2]; 
in Hi we set \vs) = [2,0,2] + [1,1,2] and in H 2 we 
set K) = [2,0,2] + [2,1,1], and (v 4 \v' 5 ) = (v' 4 \v 5 ) = 



cic\ 



"2 ■ 



Thus we obtain the following two closed chains having 
the mirror symmetry. 



Vl) 


= [4,0, 0]H 


-[3,1,0], 


v 2 ) 


= [0,0,4], 




U3) 


= [1,0,3] H 


-[0,1,3], 


v 4 ) 


= [3,0,1]H 


-[2,1,1], 


V5) 


- [2,0,2]H 


-[1,1,2]; 



(Al) 



[0,0,4] + [0,1,3], 
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K) = [4,0,0], 

K) = [3, 0,1] + [3, 1,0], 

K) = [1,0,3] + [1,1,2], 

\v' 5 ) = [2, 0,2] + [2, 1,1]. (A2) 

Even if N is very large, we can search out all these 
couples with symmetric inner product and perform the 
above procedure by a computer program. 



APPENDIX B: CALCULATION OF INVARIANT 
SUMMATION OF BASIS VECTOR INNER 
PRODUCTS 

With the [rii, ni, n^] symbol defined in section III, it is 
very easy to calculate the invariants, Sn = Yli( v i\ v 'i) = 
Si(^l^i)' f° r the pairs of closed basis vector chains, 
given arbitrary n program and m data copies in the in- 
put states. In a pair of closed chains with the mirror 
symmtry, only the terms in the form of i] (ni = n 3 ) 
contribute to this invariant sum, because |i>i)'s and |?4)'s 
carrying the terms in all other forms are orthogonal, so 
it is obtained by picking out the terms of i] satis- 
fying i < n, j < m and summing up the square of the 
coefficients absorbed in these square brackets. It can be 
done easily by a computer program for any pair of closed 
chains. 

Here we use the situation of 4 copies of program and 
3 copies of data in the input states for an example. We 
calculate the invariant sums for N = 4, 5, 6, 7 chains with 



5 elements together as shown previously. 

For N = 4 the contribution to the sum come from the 
terms [2, 0, 2] and [1, 2, 1], so the invariant sum is 

£<AM> = E^K) = %^ + ^p- = §• (Bi) 

i i 7 7 

For N — 5 the contribution to the sum come from the 
terms [2, 1, 2] and [1, 3, 1], so the invariant sum is 

/">2/"il /-il/^<3 99 

£<«} = $>M> - + = §• ( B2 ) 

i i 7 7 

For N — 6 the contribution to the sum come from the 
terms [2,2,2] and [3,0,3], so the invariant sum is 

£<*M = £teK> = ^ + ^p- = §• (B3) 

i i 7 7 

For N — 7 the contribution to the sum come from the 
terms [3, 1, 3] and [2, 3, 2], so the invariant sum is 

£<«> = £<*K> = ^ + ^ = 1- ( B4 ) 

z z 7 7 

All these pairs of chains with 5 elements have the same 
invariant sum. By induction on n and m, we see that this 
invariant sum is determined by the size of the closed basis 
vector chains. 
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